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Abstract 



["t"] ' The analytic solution is obtained describing kinetic equilibrium of 

hr] ' the /3-processes in the nucleonic plasma with relativistic pairs. The 

nucleons {n,p) are supposed to be non-relativistic and non-degenerate, 
Q^l while the electrons and positrons are ultra-relativistic due to high 

Q ■ temperature (T > 6- lO^K), or high density (p > plO^g/cm^), or both, 

H . where /i is a number of nucleons per one electron. The consideration 

^ I is simplified because of the analytic connection of the density with the 

electron chemical potential in the ultra-relativistic plasma, and Gauss 

^,.,^ . representation of Fermi functions. Electron chemical potential and 

^ I number of nucleons per one initial electron are calculated as functions 

r^ ' of p and T. 

^ ■ 1 Introduction 

o 

CN I During a collapse, leading to a formation of a neutron star, the matter passes 

states with very high temperatures and densities, and intensive birth of neu- 
trino. While in the hot core of the new born neutron star the neutrino opacity 

/S . is large enough, to establish a thermodynamic equilibrium to beta processes, 

c^ . the regions outside the neutrinosphere are almost transparent to neutrino, 

so the thermodynamic equilibrium cannot be established. It was shown in 
[B], that the characteristic time of the neutrino processes outside the neutri- 
nosphere can be much smaller than the characteristic hydrodynamic time. 
In this conditions the kinetic equilibrium is established, so that the ratio 
between neutrons and protons is determined by the conditions of the kinetic 
equilibrium, where the birth rate of protons (neutrons) in equal to its death 
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rate. The detailed investigation of the kinetic equihbrium for the pure nucle- 
onic [n—p) gas had been done numerically in [B] for a general case. Particular 
cases of the kinetic (3 equilibrium in a cold pne~ , and a hot pne^ gases have 
been considered numerically in [lOj, where the author described the kinetic 
P equilibrium approximately, in terms of the relations between the chemical 
potentials of nucleon and pairs, similar to [5j . The precision of this approach 
is rather moderate at high densities. 

Here the same problem is considered for the case of ultrarelativistic pairs 
with non-relativistic and non-degenerate nucleons. This case is considered 
without any other simplifications, and occurs to be rather simple, due to 
analytic connection of the matter density p with the chemical potential of 
the electrons /Xe, found in [2j, see also [7]. 

2 f3 reaction rates 

Let us consider kinetic equilibrium for ultrarelativistic pairs to the following 
processes 

e'+p-^n + Ue, (a) (1) 

e~^ + n —> p + Ue- (b) 

The probabilities of the reactions ([1]) for infinitely heavy nucleons are written 
as Pi 



iFti/2)n \mcC 



5 



Here (Fti/2)n ~ 1200 is the characteristic value for the neutron decay. The 
integrals I2 and J2 are defined as 



°r x^{x + xq)J{x + xqY — a^ dx 
!•? = I ■ : —. , 



'2 = / ^ 1 \ W\ ' (^) *^2) 

1 + exp(x + xq- p) 



00 



Here 



. {x + Xq + aY{x + a)\/a;2 + 2ax dx 
^ ~ J l + exp(x + a + /3) ' ^ ' 



xq = Anp/kT, Anp = rrin — nip ^ 1.29 MeV, 



a = nieC^/kT, (3 = fie/kT. (4) 

For the positrons the chemical potential /ie+ = — /^e- In the upper integral 
([3]) let us define (x = E^/kT — xq), and in the lower one (x = Ee+/kT — a), 
(i?ee+ = a/p^c^ + rriec"^), p is the electron (positron) momentum. 

In the ultrarelativistic plasma x + xq S> a in the case (a), and x ^ a in 
the case (b). Therefore the integrals ([3]) are reduced to 





/■ (x + xo)^x^ dx 



Introducing Fermi integrals 



x^(x + xo)^ dx 
1 + exp(x + Xq- l3) 

(x + xo)^x^ dx 
1 + exp(x + /3) 



(5) 



ryt "' /~f ry 

Fn{a) = / — -, (6) 



/o 1 + exp(x — a 
let us write the integrals in (jS]) as 

h = Fi{(3 - xo) + 2xoF3(/3 - xq) + x^F2(/3 - xq), (7) 

J2 = F4(-/3) + 2xoF3(-/3) + x^F2(-/3). 

In presence of ultrarelativistic pairs in thermodynamic equilibrium, there is 
an analytic expression, connecting the non-dimensional chemical potential /3 
with the matter density p, temperature T, and the number of nucleons per 
one original electron /i = "'^ "" in the nucleonic medium, written as [9l [71 [2] 



i^' + Tr'^). 




3 Kinetic beta equilibrium 

In presence of ultrarelativistic pairs the kinetic beta equilibrium is provided 
by the electron and positron captures, and the input from the neutron decay 
is negligible. In this conditions the equation of the kinetic beta equilibrium 
is written as 

T/l/(<*) T 

UpW^^^ = n^W^'\ p = l + :^^ = l + -|. (9) 

In the ultrarelativistic case thermodynamic functions as functions of p and 
T depend on the combination p/T^. Introducing a non-dimensional variable 



it follows from (|8]) and (Q the equation, determining a dependence (3{Z,xo) 
in the conditions of the kinetic beta equilibrium 

Z = (/3=^ + 7r2/3) (l + ^) , (11) 

where the integrals I2{(3,Xq) and J2{(3,Xo) are defined in (jTj). The chemical 
composition, represented by the value of /i is determined than as 

The dependences (3{p)\t and /u(p)|t in the kinetic beta equilibrium are pre- 
sented in Figs. [U-El 

The same dependences are represented in Tabs{Tll3l The results repre- 
sented in Figs. [T]-[S]and TabslT}13] coincide with numerical results from [S], 
presented in Figs 1-3 of this paper. 



4 Discussion 

Note, that the paper [S] all integrals and algebraic equations, for the compo- 
sition of the nucleonic plasma in kinetic beta equilibrium, have been solved 
numerically for a general case, while here a semi-analytical solution for the 
ultrarelativistic case is obtained. The dependence (3{p) was obtained in [10] 
using approximate approach, also for the ultrarelativistic pair conditions, the 
consideration here is exact for this case. The consideration for the kinetic 
beta equilibrium with ultrarelativistic pairs may be generalized for the mix- 
ture of nuclei in a nuclear equilibrium. In this we have several parameters 
Xzz' = ^zz'/kT. The account of the analytic connection between f3, p and p 
in the ultrarelativistic conditions, should simplify the calculations, performed 
in |1] numerically for the general case. 

The kinetic beta equilibrium is applied to the regions around the neutri- 
nosphere, in core- collapse supernovae calculations. The temperature Tj, and 
density p^ at the neutrinosphere have been calculated in [8J, giving 

T^ = 4.5 ■ 10^° - 6.3 ■ 10^° K, (13) 

p^ = 3 ■ 10" - 3 ■ 10^2 g/cml 



Table 1: The dependence of the non-dimensional chemical potential /3 = 
He/kT, and the number of nucleons per one initial electron /i = 1 + ^^, as 
a function density p, for fixed values of temperature, T = 10^° K (left), and 
T = 2- 10^" K (right), in the kinetic beta equilibrium with the ultrarelativistic 
pairs, in the nucleonic plasma, with A„p = 1.293 MeV (corresponding to 
1.5 -lO^" K). 
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Table 2: The same, as in Tab{Tl for fixed values of the temperature, T 
6 ■ 10^° K (left), and T = 2 ■ lO^^ K (right). 

~^rjK) P(g/cni^) P M T^K^ p(g/cm^) ^ J^ 
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Table 3: The same, as in TabiH for fixed values of the temperature, T 
6 ■ 10" K (left), and T = 1.5 ■ 10^^ ^ (right). 
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Figure 1: The dependence of the non-dimensional chemical potential (3 = 
fie/kT, and the number of nucleons per one initial electron yU = 1 + ^^, as a 
function density p (g/cm^), for fixed values of temperature, T = 10^° K, in 
the kinetic beta equilibrium with the ultrarelativistic pairs, in the nucleonic 
plasma, with A„p = 1.293 MeV (corresponding to 1.5 ■ 10^° K). 
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Figure 2: Same as in FiglH for for fixed values of temperature, for fixed 
values of temperature, T = 2 ■ 10^*^ K. 
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Figure 3: Same as in FiglH for for fixed values of temperature, T = 6 ■ 10^° K. 
For non-degenerate nucleons the density should not exceed 1.3 ■ 10^^ g/cm^, 
according to flTSj) . 
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Figure 4: Same as in FiglH for for fixed values of temperature, T = 2 ■ 10^^ K. 
For non-degenerate nucleons the density should not exceed 8.1 ■ 10^^ g/cm^, 
according to flTSj) . 
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Figure 5: Same as in FiglH for for fixed values of temperature, T = 6 ■ 10^^ K. 
For non-degenerate nucleons the density should not exceed 4.2 ■ 10^^ g/cm^, 
according to flTSj) . 
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Figure 6: Same as in FiglH for for fixed values of temperature, T = 1.5 • 10^^ 
K. For non-degenerate nucleons the density should not exceed 1.7-10^^ g/cm^, 
according to flTSj) . 
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The results here are apphed to the mater with non-relativistic and non- 
degenerate nucleons. The Fermi momentum of neutrons (protons) p„ (pp) in 
a fully degenerate gas is written as [2] 

/ \ 1/3 



mn,pC V6-2 ■ IQiSg/cmV 

and rrinC^ = 931.5 Me V=/c 5.4 ■ 10^^ K. It is evident from comparison with 
flT3|) . that nucleons are always nonrelativistic at the neutrinosphere. For 
nonrelativistic, nondegenerate nucleons their Fermi energy Epe should be 
less that 1.5kT. We have 

^ Uhl^ f Pn _\ ^ ^ 3 

2 ^6.2 ■ 10i5g/cmV ' ^' 2 ^ ^ 

\ 2/3 



at T = Tw = 3.6 -10'" is: 



Pr. 



v6.2-10i2g/cmV 

It follows from the comparison with flT3|) . that the temperature at the neutri- 
nosphere Tj, is always much larger than T^, so the approximation of nonrela- 
tivistic and nondegenerate nucleons is always valid near the neutrinosphere. 
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A Calculations of Fermi functions by a gen- 
eralized Gauss method 



Introducing a function 



/(^'«) = ir^^' (16) 



let us represent the Fermi function F„(a) as 

PCX) 

F^{a)= f{x,a)x"e-''dx. (17) 

Jo 

Gauss method of calculation of definite integrals suggest reducing it to an 
algebraic relation, in which it is necessary to calculate the function f{x) inside 

13 



the integral in several number m of nodes, and calculate the sum of these 
values in fixed nodes Xj, with fixed coefficients Ai, i < m. For a polynomial 
function f{x) this method gives an exact value of the integral when the power 
of the polynomial p < 2m — 1. So there is a presentation 






where the value of Xi and Ai are calculated for number of nodes between 
m = 2 and rra = 96 in [1] . To calculate improper integrals with infinite upper 
limits it is better to use a modified Gauss method, when the integrals with 
different asymptotic behavior are calculated as 

/ /(x)x"e-^dx = ^A„,/(x„,). (19) 

The values of Ani and Xni have been calculated in |3j for m = 5, n = 
0, 1, 2, 3, 4, (see also [2]). The values for n =,2 3, 4 are given in TabJH 
The Fermi functions from ([7]), with f{x,a) from (TT6l) . are represented as 

F2(a) = 0.52092 ■ /(1. 0311, a) + 1.0667 ■ /(2.8372, a) + 

+ 0.38355 ■/(5.6203, a) + (20) 

+0.028564 ■ /(9.6829, a) + 2.6271 ■ 10~^ ■ /(15.828, a), 

F-sia) = 1.2510 ■ /(1.4906, a) + 3.2386 ■ /(3.5813, a) + 

+ 1.3902 ■/(6.6270,«)+ (21) 

+0.11904 ■ /(10.944, a) + 1.2328 ■ lO'^ ■ /(17.357, a), 

Fi{a) = 4.1856 ■ /(1. 9859, a) + 12.877 ■ /(4.3417, a) + 

+ 6.3260 ■/(7.6320, a) + (22) 

+0.60475 ■ /(12.188, a) + 6.8976 ■ 10"^ ■ /(18.852, a), 
After calculation of I2 and J2 in ([7]), the values of Z, /x, T, p, are found as 



/2 ^ A Zrup (krV 
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Table 4: Roots and coefficients for calculating integrals f ll9p with m = 5, n 
2,3,4, from t2j. 



Roots Xi 








and 


n = 2 


n = 3 


n = 4 


coefBcients Ai 








Xi 


1.0311 


1.4906 


1.9859 


X2 


2.8372 


3.5813 


4.3417 


X3 


5.6203 


6.6270 


7.6320 


Xi 


9.6829 


10.944 


12.188 


X5 


15.828 


17.357 


18.852 


At 


0.52092 


1.2510 


4.1856 


A2 


1.0667 


3.2386 


12.877 


As 


0.38355 


1.3902 


6.3260 


Aa 


0.028564 


0.11904 


0.60475 


Ar, 


2.6271(-4) 


1.2328(-3) 


6.8976(-3) 
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